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ABSTRACT 

We have developed a method to calculate the fundamental parameters of the vertical 
structure of the Galaxy in the solar neighborhood from trigonometric parallaxes alone. 
The method takes into account Lutz-Kelker-type biases in a self-consistent way and 
has been applied to a sample of 0-B5 stars obtained from the Hipparcos catalog. We 
find that the Sun is located 24.2 it 1.7 (random) it 0.4 (systematic) pc above the 
galactic plane and that the disk 0-B5 stellar population is distributed with a scale 
height of 34.2 it 0.8 (random) it 2.5 (systematic) pc and an integrated surface density 
of (1.62 lb 0.04 (random) it 0.14 (systematic)) • 10~^ stars pc~^. A halo component 
is also detected in the distribution and constitutes at least 5% of the total 0-B5 
population. The 0-B5 stellar population within ^ 100 pc of the Sun has an anomalous 
spatial distribution, with a less-than-average number density. This local disturbance is 
probably associated with the expansion of Gould's belt. 

Subject headings: astrometry — solar neighborhood — Galaxy: structure — stars: 
early- type — catalogs 



1. INTRODUCTION 

Studies based on the spatial distribution of massive stars (Blaauw 1960) and interstellar hydro- 
gen (Gum et al. 1960) established some time ago that the Sun is located above the galactic plane. 
This has been confirmed using a variety of tracers (see Humphreys &: Larsen 1995 for a review) with 
somewhat incompatible values of Zq, the Sun's distance above the galactic plane, ranging from 9.5 
pc to 42 pc. Some of the discrepancies may be due to real differences in the vertical distribution of 
different stellar or gaseous components but others must be due to observational errors, since they 



^Based on data from the Hipparcos astrometry satellite. 

^The Space Telescope Science Institute is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract No. NAS5-26555. 
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appear among studies of the same type of objects. This is seen in Table 1, where the results on 
early-type stars are shown. The use of early-type stars to determine Zq has the advantage that 
their scale height with respect to the galactic plane, h, is of the same order of magnitude as Zq, 
which makes the determination simpler than for objects with larger scale heights (e.g. late-type 
stars). However, early-type stars are scarce and to obtain a statistically-reliable sample previous 
studies had to resort to objects close to the galactic plane, where extinction can introduce severe 
biases. 

The measurement of Zq using early-type stars is closely related to the vertical distribution of 
the number density of those stars, p{z). We consider here two scenarios. First, if the analyzed 
sample consists of a single population which is responsible for the gravity field in which the stars 
are moving we have a a single component, self-gravitating, isothermal disk (e.g. King 1996), which 
can be described by: 



where z is the vertical coordinate (with its origin at the Sun's position), hg is the corresponding 
scale height, and po is the density at the galactic plane. The factor of 2 is introduced in order to 
have p oc exp(— |z -|- zqI/Hs) for \z + Zq\ S> hg. A second possible scenario can be introduced by 
assuming the same conditions as before but now establishing that the gravity field is caused by a 
different population (e.g. low-mass stars or dark matter) with constant density. In that case we 
have a Gaussian vertical distribution: 



where hg is the Gaussian half-width. 

The stellar surface density, cr^t, can be obtained by integrating p{z) from — oo to -|-oo and for 
the models described by Eqs. 1 and 2 the results are 4/tsPo and V2TrhgPo, respectively. Measured 
values of h and (j+ have been included in Table 1. In each case it is indicated whether the model 
used is that of Eq. 1 or an exponential one: p{z) = po exp{—\z + ZQ\/he), in which case a* = 2hePo- 

The profiles defined by Eqs. 1 or 2 are actually not extremely different from one another and 
either one of them is expected to be a reasonably good approximation to the vertical profile of the 
disk population of early-type stars. However, we must bear in mind that there are also early-type 
stars in the halo (see, e.g. Rolleston et al. 1999 and references therein) which are expected to 
dominate the 0-B5 population at distances larger than several hundred pc from the galactic plane 
and, thus, add another component to the observed disk population. Most of these objects are 
runaway stars (Rolleston et al. 1999; Hoogerwerf et al. 2000) but some appear to have been formed 
in situ (Conlon et al. 1992). 
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The availability of the Hipparcos catalog (ESA 1997) has revolutionized the quality of astro- 
metric data and has produced for the first time good-quality trigonometric parallaxes for a large 
sample of early-type stars. However, even with these data one has to take into account a funda- 
mental problem in the use of trigonometric parallaxes for the determination of distances discovered 
by Lutz Sz Kelker (1973). Assuming a uniform spatial distribution of stars and a Gaussian distri- 
bution of observed parallaxes tTq about the true parallax tt with a width a-^^, they discovered that 
the distribution of true distances r = l/vr about the observed one Vq = I/tTq was asymmetrical and 
even divergent for r ^ oo. The effect is caused by the inverse relationship between vr and r and 
by the fact that the available volume between r and r + dr is proportional to r^. The effect is 
unimportant for parallaxes with small relative errors (Stt = (^-w/t^o ^ 0.05) and can be corrected for 
cases with e^sa 0.10 — 0.15 but for values with > 0.175 it becomes so large that an individual 
trigonometric parallax by itself provides little information about the distance to the object unless 
it is combined with additional information (e.g. the spectroscopic parallax). Lutz & Kelker (1973) 
provided a table to correct the effect of the bias for e-,^ < 0.175 (see Smith (1987); Kovalevsky 
(1998) for updated analyses) but for larger values no simple meaningful correction exists. 

A large fraction of the Hipparcos trigonometric parallaxes have > 0.175 and a non-negligible 
number have negative parallaxes. Those values cannot be discarded when statistical studies of 
distances are undertaken without introducing severe biases. One approach to the problem is that 
of Ratnatunga &; Casertano (1991), who use spectroscopic parallaxes to constrain the possible range 
of distances without discarding any stars with large values of £,r or negative ones of tTq. Another 
approach is that followed by Feast & Catchpole (1997) to evaluate the zero-point of the Cephcid 
period-luminosity relation. Those authors combine trigonometric parallaxes with up to ~ 1.0 
with the apparent magnitude and period of the stars to reduce all Cepheids to a single distance 
and produce an equivalent parallax with ~ 0.05. 

In this paper we introduce a new method to correct for Lutz-Kclkcr-type biases (section 2) 
and apply it to a sample of stars obtained from the Hipparcos catalog (section 3) in order to derive 
the parameters which determine the vertical distribution of 0-B5 stars in the solar neighborhood 
(section 4). 

2. DESCRIPTION OF THE METHOD 

The number of stars per unit distance n*^r between r and r + dr inside a solid angle J7, rii^^rdr, 
is given by: 



rii^^rdr = p{r) J7 dr, 

where p{r) is the star number density. The equivalent expression per unit parallax is: 
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If we assume that the distribution of tTq about tt is Gaussian with a width cTtt, then the 
distribution of observed parahaxes will be the convolution of the two functions: 



where ^ is a normalizing factor. Note that Eq. 4 may be singular at tt = for some choices of p{r) 
and that this may cause Eq. 5 to be non-normalizable. This is the reason why no simple correction 
can be obtained for the standard case of Lutz-Kelker bias (which has p{r) = po) when > 0.175. A 
generalization of that standard case can be achieved by allowing for non-constant density profiles, 
as previously suggested by Hanson (1979). This is motivated by the finite size of the Galaxy, 
which forces the introduction of a cutoff in any choice of p{r) based on a realistic model. Such 
a cutoff would eliminate any convergence problems in Eq. 5 and allow for finite Lutz-Kelker-type 
corrections^. 

We consider one such model based on a self-gravitating isothermal vertical distribution (Eq. 1) 
and which for early-type stars should be valid for large absolute values of the galactic latitude b, 
where the main limitation imposed on sample completeness is distance (as opposed to extinction): 



where h'g = hs\ esc 6| and zq' = Zq \ esc h\ , and the sign in the last expression is positive for 6 > and 
negative for & < 0. A least-squares fit to the observed distribution of parallaxes using Eqs. 5 and 6 
should then yield the three free parameters /i^, Zq, and (which is a function of the constants 
in the two equations). Unfortunately, two problems exist when this is attempted. First, a large 
range in h needs to be used in order to have a large enough sample, which implies that a mean 
value of I CSC 6| will have to be obtained in order to derive hg and Zq from h'g and Zq . However, for 
|6| < 25° and reasonable values of hg and zq, \ csc6| becomes so large that we are forced to go to 
distances larger than 1 kpc in order to have reasonably complete samples, and this is beyond the 
volume well-sampled by Hipparcos. Therefore, this method cannot be applied to the data currently 
available. The second problem is more subtle: Suppose we analyze all the stars contained within 
a double cone defined by |6| > bmin and that the vertical spatial distribution of early- type stars is 
given by the sum of either Eq. 1 or Eq. 2 and a halo component which is unimportant near the 



^Notc, however, that such corrections will depend on p(r) itself and that, for a non-spherically symmetric density 
distribution (as expected for a galactic disk), p(r) will be different for each direction. Thus, a tabulated version of 
Lutz-Kelker-type corrections becomes clearly impractical due to the increase in the number of degrees of freedom and 
the best strategy to estimate the distance is to compute an individual correction for each star. 
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galactic plane but dominates at high \z\. Then, since the sampled volume is proportional to z^, the 
halo component will be overrepresented in our sample and a strong bias could be introduced in our 
measurements of hs and Zq unless the proper linear combination of Pdisk and phaio is used. 

An alternative method which avoids these two problems can be designed in the following way. 
Suppose we can make an educated guess for the parameters that describe the 0-B5 disk population 
(using either Eq. 1 or Eq. 2) and that we also have a reasonably good approximation for the halo 
population (which is expected to be only a small contaminant in our sample if we avoid the bias 
mentioned in the previous paragraph). Then, for a single star with measured tTq, cttt, and b we can 
calculate the probability distribution as a function of distance as: 



where ^ is a normalizing factor and the latitude dependence is embedded in p{r) = pdisk(^') + 
Phaloif) ■ We can now project p{r) into the vertical coordinate to obtain p(z) and sum over all 
observed stars to obtain f{z), the distribution of stars as a function of the vertical coordinate. 
Doing a least-squares fit of f{z) with a vertical profile leads to a new estimate of the parameters. 
The process is then repeated until convergence is achieved. 

Several precautions have to be taken to insure the correctness of this iterative process. First, a 
maximum distance to the Sun vertical line, romaxj 

has to be specified in order to insure completeness 
and to avoid the halo-disk bias previuosly discussed. This can be achieved by truncating each 
individual f{z) at Zmax = '^maxtanft. Second, stars at very low galactic latitudes (j6| < 6min) 
have to be excluded from the sample since for lines of sight close to the galactic equator the 
spatial structure is dominated by radial (not vertical) variations and extinction. In order not to 
induce any biases due to this exclusion, the fit to the vertical profile should not include the range 
(— ■c<7maxtan6inin,^maxtan6inin)- As shown in Fig. 1, the resulting volume which will be used for 
our fit to the vertical profile is a double semi-infinite cylinder. For any given star with \b\ > 6min, 
three regions are then defined: < \z\ < tUjuax tan ^min, tUmax tan ^ \z\ ^ tUj^ax tan and 
^max tan |6| < \z\ < oo, and only the second one is used for our calculation. A star with a large value 
of will have non-negligible values of f{z) in the three regions while the probability distribution 
for one with a small value of will be concentrated in one of the three (unless its value of tTq is 
close to one of the two critical values which correspond to the two boundaries between regions). 
In any case, each of those stars will have an integrated probability in the second region between 
and 1, leading to a non- integral estimated number of stars in the double semi-infinite cylinder. 
A third precaution to be considered is the introduction of biases in the obtained parameters by 
the fitting procedure. This problem can be solved using a twofold strategy: First, biases can be 
minimized using a variable bin size determined by the criteria specified by D'Agostino &; Stephens 
(1986). Second, a population with the same parameters as those derived in the fitting process can 
be generated numerically and fitted using the same procedure. Repeating this a large number of 
times allows us to measure the numerical biases by obtaining the mean displacements between the 
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real and the derived parameters and to correct our results accordingly. 



3. 



THE SAMPLE 



We have extracted from the Hipparcos catalog the data for all stars classified therein as 0-B5 

and with |6| > &min = 5° (3 531 objects). That value of ^min was chosen in order to eliminate 
the highly obscured regions near the galactic plane while at the same time retaining most of the 
0-B5 stars in the immediate solar neighborhood. Also, for /i ~ 45 pc (a typical value chosen from 
Table 1), a vertical distance of three scale heights corresponds to ^ 1.5 kpc in m for |6| = 5°, so 
we would not expect to sample regions with large variations in the early-type star population due 
to radial galactic gradients. From that sample we have selected a subsanipic (3 382 objects) with 
good astrometric solutions by eliminating those stars with more than 5% or data points rejected 
(column 29, ESA 1997) and with a goodness-of-fit statistic g for v degrees of freedom: 



less than 3.0 (column 30)^ 

Since we expect most stars in our subsample to be located at distances > 100 pc and that for 
Hipparcos data, typically, (Ttt ~ 1 mas, Lutz-Kelker-type biases should be important. Indeed, that 
is what can be deduced from Table 2, where the distribution of values is shown. 

The other parameter which has to be selected to define the double semi-infinite cylinder which 
will be used in this work is Wj^ax- To do that, we have to consider that the Hipparcos completeness 
limit for early-type stars is determined by my = 7.9 + 1.1| sin6|, which corresponds to my = 8.0 
for |6| = 5°. A very small fraction (^ 1%) of stars below the limit is actually not included and 
the catalog is partially complete between that limit and mv ~ 12 — 13 (ESA 1997). A B5 V star 
has My —1-2; at a distance of 350 pc such a star with Ay = 1.0 mag (obtained from a very 
conservative estimate of 3.0 mag/kpc in V) would have my = 7.5, still within the completeness 
limit. If such a star was located at 135 pc (approximately 3 scale heights) from the galactic plane 
and with zu = 350 pc, its apparent magnitude would also be ~ 7.5 if we estimate Ay as ~ 0.5 
(again, a very conservative estimate at that latitude). We can then conclude that a choice of 
^i^max = 350 pc leads to a nearly complete sample for the 0-B5 disk stellar population. 



"^The selected subsample is available in electronic format at: http://www.stsci.edu/~jmaiz/data/hipparcos. 
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Table 1. Previous determinations of parameters associated with the vertical distribution of the 
early-type stellar population. The last column was calculated by multiplying the measured 
surface density by a factor derived from an unbiased 0+B sample. 



Source 


h 


Z0 


Range^ 


Model'' 


^ *, mcas 


Range"* 


^ir, equiv 




(pc) 


(pc) 






(stars pc~^) 




(stars pc~^) 


Allen (1973) 


60 




O+B 


exponential 


1.2 • 10-'^ 


O+B 


2.6 • 10"^ 


Stothers & Frogel (1974) 


46 ± 7'" 


24 ±3 


0-B5 


exponential 


5.3 • 10"^ 


0-B3.5 


0.85 • 10-3 




70 ±3*5 














Conti & Vacca (1990) 


45 ± 5 


15 ± 3 


WR 


self-gravitating 








Reed (1997, 2000) 


45 ± 20 


9.5 ± 3.5 


0-B2 


exponential 


(0.6 - 8.0) • 10""' 


0-B2 


(0.15 - 1.95) • 10-=* 



''Spectral subtypes used for the scale height and Sun's distance above the galactic plane determinations. 
''Model used for the scale height determination. 
Measured surface density. 
''Spectral subtypes used for the surface density determination. 
''Equivalent surface density for 0-B5 stars. 
'Within 200 pc of the Sun. 
BWithin 800 pc of the Sun. 



i 










^61.2 pc ^ 




700 pc 


r 


4 *■ 

r 1 



Fig. 1. — The double semi- infinite cylinder used in this paper to study the vertical distribution of 
0-B5 stars is represented by the two regions contained inside the thick lines. The values of the 
parameters used are tUmax = 350 pc and 6min = 5°. The shaded area represents the galactic disk 
with the dashed line being the galactic plane and the Sun symbol marking the position of our star. 
The hatched regions mark the area excluded to estimate the bias induced by extinction (see the 
text for details). 
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4. RESULTS 

We have carried out least-squares fits to the data described in section 3 using either a self- 
gravitating isothermal or a Gaussian disk model and applying the method described in section 
2. Since the halo population in our sample is clearly incomplete due to the Hipparcos magnitude 
limits, we used a simple reasonable model to describe it, a parabola. That choice allows us to 
discriminate between the disk and halo contributions within a few hundred parsecs of the Sun but 
can only produce a lower limit for cr^ haio- The data and fits are displayed in Figs. 2 and 3 and the 
obtained parameters are shown in Table 3. 

As previously discussed, a variable bin size was used in order to minimize numerical biases. 
For each of the two models we generated 50 populations with the same parameters derived from the 
real population and we estimated the biases by calculating the difference between the mean values 
of the derived parameters and the real ones. As expected, those biases where very small, typically 
0.5 pc for h and Zq. However, one possible bias due to extinction still remains. Suppose we are 
observing a star with b = —5° located inside the double semi- infinite cylinder but with a high value 
of £t^: This star will be counted only as a "fractional star" (i.e. its integrated probability inside the 
volume of interest will not be close to one). If the sample were complete both inside and outside 
the double semi-infinite cylinder this would not be a problem, since the lost probability would be 
compensated by that contributed by stars located outside the cylinder. However, we expect that 
extinction (and, to a lesser degree, distance) will remove from our sample some of the stars outside 
the cylinder, especially those closer to the galactic plane. To simulate this effect we consider an 
extreme case by placing a wall of infinite extinction at w = ^max which extends one scale height 
(for the self-gravitating isothermal case) or one half-width (for the Gaussian case) north and south 
of the galactic plane. The volume excluded by this wall is represented by the hatched area in Fig. 1 
for the self-gravitating isothermal case. Note that only stars in the southern galactic hemisphere 
are eliminated by this wall since the position of the Sun above the galactic plane locates any star 



Table 2. Frequency of Hipparcos relative parallax errors for the selected subsample of 3 382 
0-B5 stars used in this work. The last row refers to stars with negative parallaxes. 



£,r range 


Number 


0.00 - 0.05 


6 


0.05 - 0.15 


237 


0.15 - 0.25 


380 


0.25 - 0.50 


848 


> 0.50 


1444 


negative 


467 



-9- 



z (pc) 




Fig. 2. — Observed (histogram) and model fit frequency distributions as a function of the vertical 
coordinate for 0-B5 stars using a self-gravitating isothermal disk + parabolic halo distribution. 
The dotted lines represent each individual component and the continuous one the sum of the two. 

The displayed bin size is uniform and equal to 50 pc, but is not used for the fit (see text). The region 
immediately surrounding z = is not considered for the fit, as it corresponds to the space between 
the two semi-infinite cylinders. Note that the observed distribution shows "fractional stars" due to 
the procedure used to derive the histogram. 
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Fig. 3. — Same as Fig. 2 but using a Gaussian disk model instead of a self-gravitating isothermal 
one. 
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with 6 = 5° above the wall^. Under these conditions, we ehminated the stars in the hatched 
volume from our generated 50 populations and calculated the corresponding biases by calculating 
new least-squares fits. As expected, the new biases turned out to be more relevant than the ones 
obtained before, but even then they never induced changes larger than 15% in the values of the 
parameters. Since the real situation must be in between the two extreme cases discussed here ((a) 
all stars in the possibly excluded area partially contributing to the probability inside the cylinder 
and (b) none of them doing so), we considered the real bias to be the average of the two and we 
indicated the two extremes as a systematic error. 

The results in Table 3 indicate that both the self-gravitating isothermal and the Gaussian 
models are very good descriptions of the disk population, since g is close to in both cases and 
the fitted distribution near the galactic plane closely resembles the observed one, as shown in 
Figs. 2 and 3. The parabolic fit to the halo population is of worse quality, as expected from 
the much smaller number of stars in our sample belonging to it and from its incompleteness. 
Nevertheless, the existence of a halo population in our sample is beyond doubt, as is clearly shown 
in Figs. 2 and 3. It is also well known that this halo population extends farther than the distance- 
induced limits in our sample, even beyond a few kpc (RoUeston et al. 1999). Thus, the values of 
'7*,haio/<7*,disk which appear there have to be understood only as lower limits to the real ones. 

The values of Zq and cr^,disk obtained using the two different models are clearly compatible. 
Since we are unable to determine which of the two models provides a better fit to the disk population, 
either result can be used. Probably the ones obtained with the self-gravitating isothermal profile 
are preferred due to the better quality of the parabolic fit to the halo population, which should 
introduce less numerical noise into the result than in the Gaussian case. 

Our value for hs is somewhat lower than the ones obtained by other authors. The discrepancies 



®For the Gaussian case a small volume is also eliminated in the northern galactic hemisphere. 

Table 3. Results of the least-squares fits using either a self-gravitating isothermal (first row) or a 
Gaussian (second row) disk -|- parabolic halo distribution. The values have been bias-corrected, 
and both the random (first) and systematic (second) uncertainties are shown in each case. The h 

column refers to the scale height {hg) in the first row and to the half- width (hg) in the second 
row. No uncertainties are given for cr*,haio/<''*,disk because the fitted values are only lower- limits to 

the real distributions. 



Model 


h 


Zq 


""i.disk 


<'"*,halo/o'*,disk 




1/ 


9 




(pc) 


(pc) 


(10-3 0-B5 stars pc'^) 










1 


34.2 ±0.8 ±2.5 


24.2 ±1.7 ±0.4 


1.62 ±0.04 ±0.14 


0.058 


0.88 


41 


-0.50 


2 


62.8 ±1.7 ±4.7 


25.2 ±1.5 ±0.5 


1.67 ±0.04 ±0.24 


0.037 


0.87 


41 


-0.56 
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are only at the 2a level with the values obtained by Conti & Vacca (1990) and by Stothers & Frogel 
(1974) (within 200 pc of the Sun) and even lower in comparison with the Reed (2000) result. Allen 
(1973) does not quote an uncertainty while the value of Stothers & Frogel (1974) within 800 pc of 
the Sun is quite probably affected by heavy extinction. Some of the discrepancies may be attributed 
to contamination by halo stars (a non- negligible effect, as we have seen) and others to the use of 
a purely exponential model instead of a self-gravitating isothermal one. Our value of Zq is clearly 
compatible with the one obtained by Stothers & Frogel (1974) but is higher than the ones obtained 
by Conti & Vacca (1990) and Reed (1997), in this last case by a large amount. We believe that 
the differences may be caused by the difficulty in accounting for complex extinction variations at 
low galactic latitudes, especially when one has to rely on spectroscopic parallaxes, a problem which 
does not affect our results since they are based on unbiased trigonometric distances. This can be 
seen in the fact that the result more similar to ours is the one which is obtained from the nearest 
population. Finally, our value for o"^ is within the range obtained by Reed (2000) but is clearly 
higher than the one obtained by Stothers & Frogel (1974) and lower than the one obtained by Allen 
(1973). Once again, we believe that extinction may be responsible for some of the discrepancies 
but we must also take into account that a^, is expected to suffer significant variations depending 
on the maximum radius used, since different regions of the disk (close or far from spiral arms) are 
expected to be sampled. 

One last point to discuss is the existence of spatial variations within the region discussed in this 
paper. Stothers k. Frogel (1974) discovered that a 100 pc radius hole exists in the distribution of 
0-B5 stars at low galactic latitudes. Also, it has been known for a long time that a large fraction of 
the early-type stars in the solar neighborhood belong to Gould's belt, a structure inclined 15 — 20° 
with respect to the galactic plane and which is expanding with respect to a point located close to 
the Sun's position (see Torra et al. 2000 for a recent analysis). Is this expected minimum in the 
number density of early-type stars detected in the Hipparcos data? It indeed is. The parameters 
derived from our self-gravitating isothermal model predict that there should be 11.6 0-B5 stars 
within 66.67 pc of the Sun (i.e., with tt > 15 mas). However, only three^ 0-B5 stars are found in 
that volume: rj UMa, a Eri, and a Pav. The difference between the expected and the measured 
number appears to be too large to be explained as a statistical fluctuation, so we conclude that the 
existence of a local "hole" in the distribution of early-type stars is indeed real. 

We conclude that our method provides a solid measurement of the vertical distribution of 0-B5 
stars in the solar neighborhood by using trigonometric parallaxes alone and by taking into account 
all relevant biases. 



The three have small values of e,7 so their position can be determined with certainty to be within that volume. 
Two other stars (HIP 23342 and HIP 74368) have ttq > 15 mas but with large values of e^: and with spectroscopic 
parallaxes which clearly indicate that they are much farther away. Finally, a few stars have values of tto slightly below 
15 mas and could actually be closer to us than 66.67 pc. However, the probabilistic estimation yields a number lower 
than one, and even if there was an additional star within the given distance our arguments would remain unchanged. 
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